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RADIOMETRIC CALIBRATION AND ,? RO CESSING PROCEDURES FO R 
DEFLECTIVE BANDS ON LANDSAT-4 PROTOFLIGHT 
THEMATIC MAPPER 


ABSTRACT 


This paper provides descriptive and procedural background 
material for understanding results that are given in the 
following papers by Barker et al. appearing in these Pro- 
ceedings: '‘Prelaunch Absolute Radiometric Calibration of 

the Reflective Bands on the Landsat-4 Protoflight Thematic 
Mapper," "Characterization of Radiometric Calibration of 
Landsat-4 TM Reflective Bands," "TM Digital Image Products 
for Applications," and "Relative Radiometric Calibration of 
Landsat-4 TM Reflective Bands." 

The radiometric subsystem of NASA's Landsat-4 Thematic 
Mapper (TM) sensor is described. Special emphasis is placed 
on the internal calibrator (IC) pulse shapes and timing 
cycle. The procedures for the absolute radiometric calibra- 
tion of the TM channels with a 122-centimeter integrating 
sphere and the transfer of radiometric calibration from the 
channels to the IC are reviewed. 

The use of the IC to calibrate T) ?ta in the ground proc- 
essing system consists' of pulse .egration, pulse averag- 
ing, IC state identification, linear regression analysis, 
and histogram equalization. An overview of the SCROUNGE-era 
(before August 1983) method is presented. Procedural dif- 
ferences between SCROUNGE and the TIPS-era (after July 1983) 
and the implications of these differences are discussed. 

KEY WORDS : absolute calibration, relative calibration, band- 

to-band calibration, internal calibration system, SCROUNGE, 
TIPS, scan cycle timeline, Thematic Mapper, Landsat-4 

INTRODUCTION 

Absolute calibration is essential for a variety of scien- 
tific studies and image analysis applications. Arithmetic 
spectral transformations, such as those used to determine 
path radiance for removal or normalization of atmospheric 
effects, require absolute radiometric data. To extend sig- 
natures, based on averages or moments, beyond a scene, to 
data collected in different scenes or at different times 
under different atmospheric conditions, or even to data col- 
lected by different satellites, requires both an absolute 
measure of radiance and correction for atmospheric effects. 
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Understanding the innate bidirection reflectance character- 
istics of the target depends on exact knowledge of reflected 
radiance. To evaluate theoretical models that relate spec- 
tral radiance to ground and/or atmospheric parameteis, nu- 
merical values for radiance are necessary. 

To achieve absolute radiometric accuracy, frequent in-orbit 
recalibration o* the sensor system is necessary. The 
Landsat-4 Thematic Mapper (TM) meets this need v/ith an in- 
ternal calibration (IC) system based on incandescent lamps 
controlled by photodiode incandescent circuits. The IC sys- 
tem itself, however, is not immune to changes with time. 
Frequent recalibration by a constant external source, such 
as the Sun, would increase the absolute radiometric accu- 
racy. Use of the IC can provide relative radiometric accu- 
racy, which is essential to eliminate striping and to 
delineate boundaries. 

This paper describes the prelaunch absolute radiometric cal- 
ibration procedures and the relative radiometric calibration 
procedures in the postlaunch ground processing of the six 
reflective bands. Both SCROUNGE and TIPS (TM Image Process- 
ing System) methodologies are discussed. Lyons et al. 

(1984) present an overview of the entire SCOUNGE system, 
with special emphasis on geometric correction procedures. 
Lansing and Barker (1984) describe the calibration proce- 
dures for the thermal band (band 6) * 

THEMATIC MAPPER RADIOMETRIC SUBSYSTEM 

This section briefly reviews background essential to this 
paper and characterizes raw data obtained from the IC sys- 
tem. Engel (1980) describes the TM in detail. 

The TM observes radiance reflected from the Earth's surface 
in six bands of the following wavelengths: 0.45 to 0.52 ym, 

0.52 to 0.60 ym, 0.63 to 0.69 ym, 0.76 to 0.90 ym, 1.55 to 
1.75 ym, and 2.08 to 2.35 ym. Each band comprises 16 detec- 
tors that create 16 raster lines as the ground is scanned in 
a cross-track direction. Data are collected during the full 
cycle of the scan mirror rotation resulting in both west-to- 
east and east-to-west scans. In descending orbits (daytime), 
west-to-east scans are known as forward scans and east-to- 
west scans, as reverse- scans; in ascending orbits, the con- 
verse occurs. The 64 silicon detectors that form the first 
four bands reside in the primary focal plane. The indium 
antimonide detectors used for the remaining two reflective 
bands, bands 5 and 7, reside in a cooled focal plane. On 
each band, the eight even-numbered detectors lie in a sepa- 
rate row, staggered from the odd-numbered detectors, and 
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have separate electronic connections. Figures 1 and 2 show 
the path that light travels within the TM. 

The IC consists of three miniature tungsten filament lamps 
for the reflective bands, a blackbody for the thermal band, 
and a flex pivot-mounted resonant shutter. The IC inserts 
reference radiance levels just ahead of the spectral filters 
of each band at the prime focal plane while obscuring the 
Earth view of the detector array. To do this correctly, the 
IC shutter flag motion must be synchronized with that of the 
scan mirror assembly and the scan line corrector. (Engel 
(1980) gives a description of these subsystems.) Figure 3 
is an overview of the IC shutter flag and transfer optics. 
Figure 4 illustrates the optics for the three lamps. 

Figure 5 is a scaled schematic timeline of a TM scan cycle 
showing the time periods during which the detectors observe 
the target image and during which the shutter flag obscures 
the image during scan mirror turnaround. It should be noted 
that the image obscuration time is subdivided into a dark 
period for the collection of background data and a period 
during which the radiance from the calibration lamps is in 
view. During 3 milliseconds of the dark period, the back- 
ground radiation is adjusted to approach preset values. 
Tables 1 and 2 show all events on the TM scan timeline for a 
forward and a reverse scan. 

The lamps in the IC operate at 1900 K, with the amplitude of 
each lamp in the normal mode controlled by a silicon photo- 
diode in a feedback circuit to maintain a constant radiance 
level. (A backup mode that does not maintain the lamp radi- 
ance level is also available; with lamp aging, the backup 
mode will not provide stable lamp levels). The IC lamps, 
with energy balance filters and masks, output at levels re- 
duced by 0, 25, and 50 percent (see Figure 4). A bundle of 
six fiber optic threads delivers the combined lamp output to 
the shutter flag. This method of transmitting radiation to 
the moving calibration shutter allows the lamps to provide 
radiation independently and to contribute proportionately 
toward illuminating all detectors. When in the automatic 
sequencing mode, the three lamps are cycled automatically to 
provide eight calibration levels (seven plus dark) to all 
detectors of bands 1 through 5 and 7. Each level is on for 
approximately 40 scans. In this paper, lamp configuration 
is indicated by three binary digits (that is, 100 means 
lamp A is on and B and C are off, Oil means only B and C are 
on , and so on) . 

The shape of the pulse output by the internal calibration 
system depends on the physical location of the detector with 
respect to the calibration light source. Detectors are 
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FIGURE 5 

DIGITAL MINOR FRAME LOCATIONS 
FOR SUCCESSIVE TM SCANS 
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FOOTNOTES FOR TABLE 1 

MF LOCATIONS FOR ODD-NUMBERED (FORWARD) SCANS 


(1) Collection of digital data for radiometric calibration of TM during- 
the Scrounge-Era preprocessing by ADDS, f|4mf (minor frames) of 1C 
(Internal Calibrator) data are temporarily put in a buffer of a Macro Array 
Processor (MAPI). These calibration data are collected starting at minor 
frame 6407, which is close to the end of the scene video data and the 
beginning of shutter obscuration. 200 of the 824mf are sent on to a Vax 
computer for use in radiometric calibration. 

For odd-numbered scans (forward sweeps of the TM scan mirror), these 200mf 
from the IC shutter are collected from three separate regions (windows). 
Labels and descriptions for the three collect windows from each forward 
scan are given below: 

BF-BDC This is a 24mf region of dark level background (BKG) taken 
on a forward scan before DC restoration begins. Prior to 
May 20, 1983, there was a single 52mf ADDS BKG collect 
window which started before DC restoration, at MF 6643. 

BF-ADC 28mf BKG taken after DC forward restoration ends. 

CALFWD This is a 148mf forward scan calibration (CAL) region 
containing the TM responses to light from the current 
configuration of the three IC lamps. Prior to December 22, 
1982, the 148mf CAL collect window started at MF 7009. It 
then, began at MF 7029 until May 20, 1983, when it was 
changed to MF 7037. 

(2) MF (Minor Frame) locations were arrived at one of four ways: 

DEF By Definition 

NOM From nominal value of another variable 

OBS Observed from in-orbit digital data 

CALC Calculated from values of other MF locations and defined, 
nominal or observed relative differences 

(3) MSEC (Milliseconds) locations were calculated from MF locations 
assuming a nominal 9.61) microseconds per minor frame. 

(4) TRAPP is a TM Radiometric and Algorithmic Performance software program 
run on pre and postlaunch tapes at the LAS Facility. 
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Table 2 

Timing Locations on Hain Shutter during Even Numbered (Reverse) Scans 
of Internal Calibrator (IC) on Landsat-4 Thematic Mapper 
and Scrounge-ERA (1982-1383) Ground Processing IC Collect Windows 
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FOOTNOTES FOR TABLE 2 

MF LOCATIONS FOR EVEN-NUMBERED (REVERSE) SCANS 


(1) Collection of digital data for radiometric calibration of TM during 
the Scrounge-Era preprocessing by ADDS, 824 mf (minor frames) of IC 
(Internal Calibrator) data are temporarily put in a buffer of a Macro Array 
Processor (MAPI). These calibration data are collected starting at minor 
frame 6407, which is close to the end of the scene video data and the 
beginning of shutter obscuration. 200 of the 824mf are sent on to a Vax 
computer for use in radiometric calibration. 

For even-numbered scans (reverse sweeps of the TM scan mirror), these 200 
mf from the IC shutter are collected from three separate regions (windows). 
Labels and descriptions for the three collect windows from each reverse 
scan are given below; 

CALREV This is a 148mf reverse scan calibration (CAL) region 
containing the 7M responses to light from the current 
configuration of the three IC lamps. Prior to December 22, 
1982, the 148mf CAL collect window started at MF 6626. It 
then began at MF 651,7 intil May 20, 1983, when it was changed 
to MF 6513, 

BR-BDC This is a 28mf region of dark level background (BKG) taken on 
a reverse scan before DC restoration begins. 

BR-ADC 24mf BKG taken after DC reverse restoration ends. Tnis is 
the same shutter location as the 24mf forward BKG. Prior to 
May 20, 1983, there was a single 52mf ADDS BKG collect window 
which started after DC restoration, at MF 7089. 

(2) MF (Minor Frame) locations were arrived at one of four ways: 

* 

DEF By definition 

NOM From nominal value of another variable 

OBS Observed from in-orbit digital data 

CALC Calculated from values of other MF locations and defined, 
nominal or observed relative differences. 

(3) MSEC (milliseconds) locations were calculated from MF locations 
assuming a nominal 9.611 microseconds per minor frame. 

i 

(4) ' TRAPP is a TM Radiometric and Algorithmic Performance software program 
run on pre and postlaunch tapes at the LAS Facility. 
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stored by band number in even-numbered and odd-numbered 
rows; thus even (or odd) numbered channels of a given band 
have the same pulse shape. Sample pulses are provided in 
Appendix A (Figures A-l through A-18) ; Figure 6 shows a 
typical pulse. Pulse shape is independent of scan direc- 
tion. If the pulse from a reverse scan is overlaid on a 
forward scan pulse, the difference in counts for correspond- 
ing positions is at most two digital numbers (DN) (Figures 7 
and 3) . 

Within a given band, the location of the pulse center prog- 
resses to higher numbered minor frames as the channel number 
decreases. Bands 1 and 4 have sharp asymmetric peaks; 
bands 2 and 3 ha*'* rounder, less asymmetric peaks. Bands 1 
and 2 have peaks to the right of the pulse center in the 
forward scan; bands 3 and <1 have peaks to the left of the 
pulse center. In bands 1 and 2, the odd channels have 
larger pulses, about 6 counts higher in band 1 and 20 counts 
higher in band 2. The even channels in bands 3 and 4 have 
slightly larger pulses than the odd channels, about five 
counts in band 3 and two counts in band 4. Band 5 even 
pulses are 38 counts higher than odd pulses. For both even 
and odd channels, the pulse is asymmetric, peaking to the 
right of the center, but the odd-channel pulses are much 
rounder. Although even-channel pulses in bands 6 and 7 are 
only <• ight counts greater than odd-channel pulses, the pulse 
shap>. £ #re qualitatively different in the two cases. The 
ev&i channels have nearly symmetrical rounded pulses, peak- 
ing slightly to the right of center. The odd channels have 
a flat plateau to the left of center, and then a peak to the 
right. 

Figure 9 shows pulse averages, as determined by the Hughes 
method discussed below, versus scan number for a post-launch 
scene. It should be noted that transitions involving the 
turning on of a lamp result in approximately 15-percent 
overshoot. In mid-IR bands, bands 5 and 7, blackbody heat 
radiation from the newly turned off lamp contributes to the 
pulse height for approximately eight scans (four forward and 
four reverse) . 

ABSOLUTE RADIOMETRIC CALIBRATION PROCEDURE 

Prelaunch absolute radiometric calibration of the TM sensor 
channels and the IC is critical because the TM does not 
observe the Sun or other extraterrestrial bodies of known 
radiance while in orbit. This section presents an overview 
of prelaunch absolute calibration procedures used on the 
Landsat-4 TM. A detailed discussion of the results of ab- 
solute calibration is presented in Barker et al. (1984a). 
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RADIANCE. Q (DIGITAL COUNTS) RADIANCE. Q (DIGITAL COUNT® 


ORIGINAL PAGE! !$• 
OF POOR QUALITY 


FIGURE 6 

RADIOMETRIC CALIBRATION - TM LANDSAT-4 
TM1 d FORWARD SCAN, MARCH 9, 1382, VACUUM PRELAUNCH 



FIGURE 7 

COMPARISON OF CALIBRATION PULSES BETWEEN FORWARD 
AND REVERSE SCANS PRELAUNCH DATA 



12 









The prelaunch radiometric calibration procedures are de- 
signed to ensure that the TM meets its performance require- 
ments. The radiometric specifications can be divided into 
three groups: dynamic range, sensitivity, and accuracy. 

The TM is designed for land observations? the reflective 
bands are therefore designed for maximum dynamic range for 
the expected land radiance levels (Fraser, 1975; Duck, 

1977) . This range can be tested using a well-calibrated 
external source. The radiometric sensitivity is expressed 
in terms of s ignal-to-noise ratio (SNR) for reflective 
bands. Radiometric accuracy is one of the most important 
issues to be addressed. The relative accuracy requirements 
are quite stringent: about 1 quantum level for between- 

detector accuracy and 2 percent for between-band accuracy. 
The overall absolute accuracy specified for the TM is 
10 percent. Obtaining absolute calibration accuracy for 
spaceborne instruments has traditionally been a very diffi- 
cult tank. 

The calibration procedure used was developed by Hughes 
Aircraft Company (HAG) to test the performance of the TM 
channels; to determine the end-to-end radiometric transfer 
function with 10-percent accuracy; and to verify that the 
relative radiometric accuracy is better than 2 percent, the 
within-band accuracy is 1 DN, the SNR is within specifica- 
tion, and the IC is accurate and has the correct range. 

CALIBRATION EQUIPMENT 


The usual approach to the prelaunch calibration of visible 
and near-infrared sensors is by comparing them to a sec- 
ondary source that has been calibrated with a National 
Bureau of Standards (NBS) standard lamp. The secondary 
standard used with the TM is a 122-centimeter (diameter) 
Integrating Sphere (IS(122)) maintained by the National 
Aeronautics and Space Administration/Goddard Space Flight 
Center (NASA/GSFC) and HAC Santa Barbara Research Center 
( SBRC ) personnel. SBRC has maintained the IS (122) and has 
periodically recalibrated it by comparison with a secondary 
standard lamp using a monochromator. These tests were per- 
formed in June 1979 and February 1982 at NASA/GSFC and in 
June 1980 and May 1982 at SBRC and are discussed by walker 
(1982a, 1982b). The secondary standard is related to an NBS 
standard by a similar calibration procedure. 

The IS (122) is a hollow sphere with an inside diameter of 
122 centimeters. It consists of two aluminum hemispheres 
mounted in a framework that allows easy access to the in- 
terior for bulb and surface maintenance. A 46-centimeter- 
diameter aperture allows reflected light froa the 
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12 quartz-halogen standard lamps out of the interior to cal- 
ibrate the sensor. The inside surface of the sphere is 
coated with a white diffusely reflecting barium sulfate- 
based paint (McCullough et al., 1969). 

The lamps in the IS(122) have different wattages. There are 
six 200-watt lamps, two 100-watt lamps, and four 25-watt 
lamps, for a maximum of 1500 watts. The lamps are operated 
at derated wattages to prolong their useful lifetimes. 

During the calibration tests, the lamps are all lighted 
initially for each test sequence, as shown in Table 3. The 
lamps are then sequenced off to reach the next lower radi- 
ance level. Twenty levels are used in all. The responses 
of all 96 detectors in the reflective bands (bands 1 through 
5 and 7) are noted at each level, and the detector gains, 
offsets, and SNRs are calculated. 

The IS (122) is directly calibrated against a secondary 
standard lamp that has in turn been calibrated by a standard 
lamp of similar design and is directly related to an NES 
standard. The standard and the secondary lamps are both 
quartz-halogen 1000-watt lamps. The overall accuracy of 
these lamps and the complete calibration sequence are dis- 
cussed in Barker et al. (1984a) . 

C ALIBRATION PROCEDURES FOR THE TM WITH THE IS (122) 

Calibration of the TM with a known source, the IS(122), con- 
sists of comparing the TM output with the series of known 
radiances input to the TM from the source. This procedure 
is designed to determine the overall response of the TM 
under known conditions. 

Each calibration of the IS (122) results in a series of spec- 
tral radiance values at intervals of 50 nanometers in wave- 
length from 400 to 2500 nanometers for each sphere-lamp 
configuration (combination of particular lamps on) (Walker, 
1982b). This spectral radiance output from the IS (122) is 
then combined with the measured wavelength response of the 
TM (Markham and Barker, 1984) in each band to obtain the 
spectral radiance that each TM detector would see in a par- 
ticular band and for a particular IS(122) lamp configuration. 

These radiance levels by band are the known input to the TM 
for calibration. During each test, the digital output re- 
sponse of the TM to each input level is collected, averaged, 
and then compared to the known radiance input to obtain the 
calibration parameters (gain and offset) . Figure 10 illus- 
trates a linear fit of digital counts to radiance to deter- 
mine channel gain (G(C)) and offset (0(C)). 
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INTEGRATING SPHERE LAMP CONFIGURATIONS BY RADIANCE LEVEL 
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The SNR at each^level (k) is also computed from the averaged 
digital count (P^) and, the standard deviation ( 0 |<) : 


SNR 



A least-square linear regression is used to obtain SNR as a 
function of P. 

The input radiance to the TM is computed by combining the 
output of the IS (122) for a given configuration of lamps 
(radiance level) with the measured response of the TM in 
each band. This is given by 


L x< B > 



L 


j 

sphere 


(X) 


R(X) dX 



2 

. R (X ) dX . 


where 

L^(B) = spectral radiance in a particular TM band 
« • 

L^ph ere (X) = output of the sphere at radiance level j 

R (X ) ~ spectral response of the TM band as a 
function of wavelength 

The integration is done from wavelength X]_ to X 2 for 
a particular TM band. The critical parameters in the 
calibration are alignment and separation between the IS(122) 
and the TM, the TM scan mirror operating mode (locked or 
scanning) , and the sphere output to the TM in each band at 
each radiance level (L^(B) above). The alignment of the TM 

with the IS (122) is shown schematically in Figure 11. The 
distance marked A in the figure has varied widely in the 
different tests performed to calibrate the TM. Table 4 gives 
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the dates, the separation distances (A) between the TM and 
the IS (122), and a letter designation for each test. 

Table 5 shows the variability in the measured channel gain 
in each test. The alignment of the TM with the IS (122) is 
particularly important at the larger separations if the scan 
mirror is locked in place. The IS ( 122 ) aperture is 
46 centimeters, which is illuminating the 41-centimeter TM 
aperture; at a separation of 5.2 meters, the spread of the 
TM spectral bands along the scan is 3 centimeters at the 
prime focal plane. This leaves a maximum travel of 2 centi- 
meters during which all detectors of each band (1 through 5 
and 7) may simultaneously see the same part of the IS (122). 
If the scan mirror is locked in place, the alignment thus 
becomes critical. This is not the case if the scan mirror 
is allowed to move. Test procedures originally required the 
scan mirror to be locked at midscan; these were changed 
after the June 1981 test. Both locked and scanning modes 
were used in March 1982. To ensure that the peak response 
is actually measured, the test procedures with the scan mir- 
ror operating may also require that the raw output from the 
sphere be searched for the largest response and then fitted 
by a second-order polynomial. This has not been a part of 
the calibration procedure, however, and may have to - be in- 
cluded to reduce the error if it is found that the output of 
the IS ( 122 ) as seen by the TM is not sufficiently flat. The 
test procedure may thus be summarized as follows: 

1. Select an external radiance source as standard 
(IS (122) ) 

2. Determine the characteristics and output of the 
external standard 

3. Adjust trim resistors for each detector in the TM 
to obtain the correct dynamic range 

4. Measure the TM performance at various external 
radiance levels and cooled focal plane temperatures 

5. Determine the gain, offset, dynamic range, and SNR 
of each detector 

6. Transfer the calibration from the detectors to the 
IC 

7. Use the TM to calibrate other external calibrators 
for use when the IS(122) is unavailable (as in 
vacuum chamber testing) 

8. Conduct spectral matching tests for band-to-band 
comparisons 
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TABLE 5 

SUMMARY OF AVERAGE GAIN CHANGES IN PERCENT BY TEST 
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CALIBRATION TRANSFER TO THE IC 


The IC is used mainly for measuring relative performance in 
flight. External calibrators can be used only on the 
ground. To maintain good radiometric accuracy, it is impor- 
tant to carefully calibrate the IC using a known standard. 
Furthermore, the behavior of the TM as a function of temper- 
ature and age should be well understood so that IC measure- 
ments can be corrected, because there is no way to obtain 
absolute calibration in flight. 

The function of the calibration transfer is to provide ref- 
erence radiance levels for the IC in different operating 
modes. These reference data are then used to calibrate in- 
orbit scene data and to monitor tne health of the detectors 
during postlaunch operation. The procedures used to process 
the raw transfer data should be the same as those proposed 
for operational postlaunch ground processing. To transfer 
the absolute IS (122) calibration from the TM detectors to 
the IC, the raw IC calibration pulse must be extracted from 
the telemetry stream; a suitable pulse integration method 
must be applied to all detector data; and all integrated 
values for a particular lamp configuration must be averaged, 
excluding certain scans near lamp state changes. The actual 
algorithms used are the same as the SCROUNGE ground radio- 
metric processing procedures discussed in the next section. 
The lamp state values may then be converted from averaged 
digital counts (P) (output from the TM) to effective spec- 
tral radiance using the calibration parameters (G and 0) de- 
termined from the IS(122) calibration (Barker, 1984), 

Regardless of the IC operational mode, all eight lamp states 
for all detectors must be processed through the eight steps 
mentioned previously. 

GROUND RADIOMETRIC PROCESSING PROCEDURES 

From launch through July 1983, TM image data were processed 
using the IC with the SCROUNGE processing algorithms. In 
July 1983, an early version of TIPS became available. The 
important features of the SCROUNGE algorithms will be de- 
scribed here, and the resulting products will be mentioned. 
TIPS will be discussed only in terms of its differences from 
SCROUNGE. 

The use of the IC to radiometr ically calibrate image data 
from the reflective bands involves six steps: 

1. Extracting calibration minor frames 

2. Integrating a pulse to obtain a pulse height 
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3. Averaging pulse heights belonging to a state and 
identifying the IC state 

4. Obtaining detector gains and offsets by regression 

analysis of observed pulse averages versus nominal 
values . 

5. Adjusting gains and offsets of channels within a 
band to achieve a common radiance range 

6. Balancing relative grey levels using histogram 
normalization (optional) 

During the SCROUNGE era, only the minor frames (mf) falling 
within a 148-mf window are examined for pulse extraction. 

If part of the pulse falls outside the window, the inte- 
grated pulse value is reduced. The pulse extraction window 
is located at a fixed distance from the start of the line. 
The distance is, the same for all channels but different for 
forward and reverse scans, A detailed history of the loca- 
tion of calibration windows is given in Barker et al. 
(1984b), SCROUNGE integrates pulses by using the Hughes 
algorithm (Figure 12) , which consists of the following steps 

1. Locate pulse maximum 

2. Move outward from the maximum to locate pulse edges 
that are 40 percent of the maximum value 

# 

3. Define pulse center as the midpoint of the line 
segment defined by the pulse edges 

4. Average 31 points (pulse integration width) about 
the pulse center to obtain the pulse height 

Pulse averaging involves identifying an IC statf transition, 
determining which scans between transitions to include in 
the average, averaging the pulses from selected scans to 
obtain the averaged pulses P(A,C), and associating 
P(H r C) with an IC state. In SCROUNGE, pulse heights from 
forward and reverse scans are averaged together, A 
transition occurs when the pulse height of a scan differs 
from the pulse height of the previous scan by more than 
three digital counts. The first 25 scans following a 
transition are ignored (Figure 9) ; the remaining scans, 
preceding the next transition, are averaged to determine Q. 
The average pulse is associated with theIC state having the 
closest nominal averaged pulse value. Data from each of the 
96 reflective channels are processed separately. When 






51 channels agree on an IC state, the XC state is considered 
identified and is assigned to all channels. 

Gain (G (G) ) and offset (0(C)) are determined for each 
channel by a linear unweighted regression of counts versus 
nominal counts using all eiyht IC states. A radiance range 
in digital counts common to all the channels within the band 
(RMAX-RMIN) is determined for each band. Gains and offsets 
are then normalized so that values obtained for each channel 
lie within the common range. 

First RMAX and RMIN are used to determine postcalibration 
gain (G°(B)) and offset (0°(B))s 


G° (B) 


QMAX 

RMAX-RMIN 


n °, m _ QMIN - QMAX*RMIN 
' ' (RMAX-RMIN) 


where QMAX is the full-scale digital value (255) and QMIN is 
the minimum digital value (0). Normalized gain NG IC (C)and 


normalized offset (NOj ( -,(C) ) for each channel (C) in the band 
are then determined as follows: 


NG JC (C) - G(C)/G°(B) 

NOjc(C) = 0(C) - NGxc(C) * 0°(B) 


As a result of this approach, adjusted count values for par- 
ticular channels that lie outside the common range are 
simply set to lie at the maximum or minimum of the range as 
appropriate. 

Histogram normalization involves the following steps for 
each band: 

1. Histograms are obtained for each channel. 

2. The histograms for all the channels in the band are 
combined to form an average histogram. 
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3 . 


A new gain (N'G) and a new offset (N'O) for each 
channel are derived by adjusting NG IC (C) and NO IC (C) 

so that the channel histogram will have the same 
average (p) and standard deviation (a) as the 
band average histogram: 


N'G = 


NG IC (C) 


and 


a * p. 

N'O = NO IC (C) - N'G * £ 


where Pi and <jj[ are the respective mean and 
standard deviation for the histogram for channel i 
before normalization. 

4. The new gains and offsets are then used to produce 
a radiometric lookup- table (RLUT) for each channel: 


RLUT... 



- N'O 
n t g 


where RLUT^j is the value in the lookup table 
corresponding to an initial value Pj in channel C. 

5. The image is radiometr ically corrected using the 
RLUT. 

Uncalibrated image data are available on a computer- 
compatible tape (CCT) known as the CCT-BT tape. This tape 
also contains the initial gains and offsets (prior to range 
adjustment and histogram normalization) for each channel. 
CCT-BT tapes generated after April 1983 also contain the 
histogram-normal: "ed gains and offsets. The CCT-AT tapes 
contain the radiometr ically corrected image. A final prod- 
uct containing radiometr ically and geometrically corrected 
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resampled data is also available. These digital products 
are described in detail in Barker and Gunther (1984) . 

TIPS is a flexible system that provides three options: the 

use of nominal gains and offsets, the use of nominal gains 
and scene-specific background for offsets, and the use of 
the IC to determine gains and offsets. In each case, 'the 
use of histogram normalization is optional. Differences 
between the SCROUNGE and TIPS processing algorithms are 
highlighted below. 

In general, the TIPS algorithms produce more stable results 
and handle noise in the calibration data better. Specifi- 
cally, TIPS will examine the entire shutter flag part of the 
scan to extract the calibration pulse, thereby eliminating 
the risk of clipped pulses. The TIPS pulse integration 
method, THRESH, uses a wider pulse integration width of 
65 minor frames. The pulse averaging algorithm uses infor- 
mation from several scans (currently 10) to locate IC state 
transitions. Out-of-range scans (deviating by a predefined 
amount from the average of the stable part of the current 
data) as well as scans on both sides of a transition are not 
included in the averaged pulse value P. In the regression 
analysis, background data from each scan as well as the 000 
IC state value are used. Each pulse average used in the 
regression is weighted based on the number of scans contrib- 
uting to the pulse average and the standard deviation of the 
average. The one aspect of SCROUNGE that is probably safer 
than TIPS is in IC state identification. Unlike SCROUNGE, 
TIPS uses only one or a small number of channels to identify 
the IC state. 


SUMMARY 


The TM radiometric subsystem observes radiance reflected 
from the Earth in six reflective bands: 0.45 to 52 pm, 

0.52 to 0.60 pm, 0.63 to 0.69 pm, 0.76 to 0.90 pm, 

1.55 to 1.75 pm, and 2.08 to 2.35 pm and one thermal 
band at 10.4 to 11.6 pm. Each band comprises 16 detectors 
that create 16 raster lines as the ground is scanned in both 
forward and reverse directions along the crosstrack. The IC 
consists of three miniature tungsten filament lamps for the 
reflective bands, a blackbody for the thermal band, and „ 
flex pivot-mounted resonant shutter. The IC inserts refer- 
ence radiance levels just ahead of the spectral filters of 
each band at the prime focal plane while obscuring the Earth 
from the detectors. 
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Prior to launch, the TM channels were absolutely calibrated 
with a secondary standard 122-centimeter integrating 
sphere. The objectives of the calibration procedure were to 
determine absolute channel gain and to characterize relative 
radiometric accuracy, within-band accuracy, SNR, and dynamic 
range. The TM channels were used to calibrate the IC. 

t 

The use of the IC to calibrate TM data in the ground proc- 
essing system consists of pulse extraction, pulse integra- 
tion, pulse averaging, IC state identification, linear 
regression analysis, adjustment of the detectors within a 
band to a common range, and histogram equalization. The 
important features of the SCROUNGE radiometric calibration 
system, which was operative through July 1983, are as 
follows : 

• Pulse integration width is 31 minor frames wide. 

• Forward and reverse scans are processed together. 

• The first 25 scans following an IC state transition 
are excluded from the pulse average. No scans are 
skipped before a lamp transition. 


• IC state identification is made through comparison 
with nominal values. All channels are checked with 
the results determined by agreement of the majority. 


• Linear regression analysis is unweighted and in- 
cludes all eight IC states. 


• In histogram equalization, the mean and standard 
deviation of each channel are adjusted to conform 
with the band average. 


The TIPS radiometric calibration process, which replaced 
SCROUNGE in July'1983, differs from SCROUNGE as follows: 


• Pulse integration width is 65 minor frames. 


• In pulse averaging, scans on both sides of an IC 
state transition are skipped. 

• Only a small preselected group of channels are 
checked to identify the IC state. 

• Linear regression is weighted. Background data 
from each scan are used as well as the eight IC 
states. Ir' band 4, the 111 state, which is satu- 
rated, is n it used in the calibration. 
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RECOMMENDATIONS 


Brief recommendations are given below; expanded recom- 
mendations are given in Barker (1984) „ An improvement in 
both relative and absolute accuracy of the radiometric 
calibration can be achieved by modifying the calibration 
procedure. Of course/ any change in the postlaunch pro- 
cedure would require a reprocessing of the prelaunch 
absolute calibration data using the modified algorithms. 

Prom the results presented in Barker et al. (1984b) , it is 
clear that forward and reverse scan data should be processed 
separately. h wider pulse integration width would reduce 
the noise in the averaged pulse values. A wider pulse ex- 
traction window or a pulse extraction window that shifted 
with channel number would eliminate pulse truncation and 
permit the use of a wider integration width. The use of 
background data from each scan line for calibration of that 
scan would probably be preferable to applying the average 
background observed during the lamps-off (000) state to all 
the scans in a scene. When applying background data to an 
image line/ however/ care must be taken to use background 
data that have not been altered by dark current (DC) restore 
subsequent to the time the image was taken. Barker et al./ 
1984/ presents additional recommendations for postlaunch 
procedures; Barker et al., 1984a, recommends improvements 
for the absolute calibration procedures. 
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APPENDIX - CALIBRATION PULSES 


Figures A-l through A-18 are calibration pulses taken from 
the March 9, 1982, prelaunch vacuum data. Samples are in- 
cluded for forward and reverse scans, even and odd channels. 
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